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EXPERIMENTAL INVESTIGATION OF A TURBULENT 


BOUNDARY LAYER AT MACH 6, HIGH REYNOLDS NUMBERS, 

AND ZERO HEAT TRANSFER 

By Jerry B. Adcock, John B. Peterson, Jr., 
and Donald I. McRee 
Langley Research Center 


SUMMARY 

An experimental investigation of the turbulent boundary layer on a hollow 
cylinder was made at a Mach number of 6 and with adiabatic wall temperature at 
Reynolds numbers based on the distance from the leading edge from 5 X 10^ to 
33 X lo6. Data obtained include turbulent recovery factors, a total tempera- 
ture profile, velocity profiles, and average skin-friction coefficients by the 
momentum method. 

The turbulent recovery factor was found to be approximately 0.88; a slight 
decrease of recovery factor with increase in Reynolds number was obtained. The 
velocity profiles, when compared with power profiles, were found to have an 
increasing exponent with increasing Reynolds number. The measured temperatures 
in the boundary layer were compared with two theoretical temperature distribu- 
tions, the Crocco laminar distributions used by Van Driest (NACA TN 2597)/ a 
quadratic distribution used by Donaldson (NACA RM L52H04) , and with the assump- 
tion of a constant total temperature through the boundary layer. The assump- 
tion of a constant total temperature was found to give the most accurate value 
of momentum thickness. The skin-friction results were compared with the 
theories of Van Driest, Wilson and Van Driest, Cope and Monaghan, with the 
Sommer and Short T* method, and also with other experimental data. 

INTRODUCTION 




Since turbulent boundary layers will exist on large portions of the surface 
of high-speed vehicles, it is essential to know their characteristics. Turbu- 
lent boundary layers, however, are not as amenable to theory as laminar boundary 
layers; therefore, experimental determination of their characteristics must be 
made in order to check the validity of the various theories. In turbulent flow 
there is wide variation in the many theories. This variation can be seen in 
figure 1 of reference 1, in which the large difference in values predicted by 
the various skin-friction theories is shown. Many correlations of the skin- 
friction theories and the available experimental data have been made. One of 


the latest correlations can he found in reference 2. It should be noted, 
however, that most of the data were obtained at the lower Mach numbers and' 
Reynolds numbers and that the differences in the skin-friction theories become 
greater at the higher Mach numbers. 

The purpose of this investigation was to provide turbulent -boundary-layer 
data at a hypersonic Mach number of 6 and high Reynolds numbers based on the 
distance from the leading edge (5 X 10& to 33 X 10^ ). Turbulent-boundary-layer 
velocity profiles were obtained from pitot-pressure surveys through the boundary 
layer on a hollow cylinder. These velocity profiles are compared with theoret- 
ical power profiles. Average skin-friction coefficients were obtained by use 
of the momentum method. In this method the momentum thickness is calculated 
from a pitot-pressure survey and from a temperature distribution through the 
boundary layer. Normally, at low Mach numbers the assumption of a constant 
total temperature through the boundary layer is made. In order that a better 
insight into the temperature distribution which should be used at a Mach number 
of 6 be gained, a total-temperature survey was made. 

In addition to the skin-friction coefficients and velocity profiles calcu- 
lated from the temperature and pitot -pres sure survey data, turbulent recovery 
factors were also obtained from wall-temperature measurements. 

The present skin-friction coefficients, along with other authors' experi- 
mental skin-friction data, are compared with four of the well-known theories: 
Sommer and Short (ref. 3) > Van Driest (ref. 4), Wilson (ref. 5 or 6) and 
Van Driest (ref. 7) ) Cope (ref. 8), and Monaghan and Johnson (ref. 9) • 


SYMBOLS 


A,B , C ,D,K, 
a,b,c 


constants 


Cp average skin-friction coefficient 

Cf local skin-friction coefficient 

Cp specific heat at constant pressure 

i limits of summation 

M Mach number 

Mg Mach number in constant-pressure plane (average Mach number at 

edge of boundary layer) 


n 


exponent in power-profile relation, — 

u 5 
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pressure 

pitot pressure 

Reynolds number per foot 

Reynolds number based on distance from leading edge of model 

Reynolds number based on distance from leading edge to virtual 
origin of turbulent boundary layer 

Reynolds number based on distance from virtual origin of turbulent 
boundary layer 

Reynolds number based on boundary-layer momentum thickness 
Reynolds number based on 0 

radius of cylinder, 3 inches 

temperature 

total temperature 

velocity at outer edge of boundary layer 
velocity in boundary layer 

distance from virtual origin of turbulent boundary layer 
distance from leading edge of model 

distance, or height, normal to model, measured from surface 

angle of attack 

ratio of specific heats 

boundary- layer thickness 

boundary-layer thickness based on pitot-pressure profile 
boundary- layer displacement thickness (see eq. (9)) 
boundary-layer shape factor 
recovery factor 

momentum thickness of boundary layer (see eq. (8)) 
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0 

0 H 

e^e 11 

p 


momentum thickness based on Mach number in plane of constant 
pressure (see eq. (10)) 

boundary-layer total-enthalpy thickness (see eq. ( 3 )) 

coefficients used in equation (6) 
gas density 
angle of yaw 


Subscripts: 

aw 

i 

max 

P 

t 

w 

5 


adiabatic wall 
incompressible 
maximum 

measured by boundary-layer probe 
total or stagnation conditions 
wall or wall conditions 
conditions at edge of boundary layer 


APPARATUS 


Wind Tunnel 

The tests were conducted in the Langley 20-inch Mach 6 tunnel. This tunnel 
is of the blowdown type; however, the capacity of the air-supply pumps was large 
enough to allow test runs as long as 45 minutes in this investigation. The 
tunnel has a rectangular test section 20 inches by 20.5 inches. Maximum oper- 
ating stagnation pressure is about 525 psia and maximum stagnation temperature 
is about 600° F. A schematic drawing and further information can be found in 
reference 10. 


Model 

The hollow-cylinder model used in the test is shown schematically in fig- 
ure 1. The model outside diameter is 6 inches for a length of 42 inches and 
then flares out at an angle of 20° to a diameter of 8 inches. This larger 
diameter was needed to facilitate the circulation of a coolant between the 
inside and outside walls of the model for further cooled-wall boundary- layer 
investigation. The leading edge was beveled inward at an angle of 15° until 
the inside diameter of the model was 5 inches. 
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Located at 1.15 inches behind the leading edge was a transition trip used 
to cause early transition. This transition trip was made of 0.075-inch- 
diameter steel rods mounted normal to the surface by pressing them into drilled 
holes. The trip was made with the rods spaced every 0.25 inch around the cir- 
cumference of the model with an original height of 0.075 inch. A preliminary 
study using fluorescent-oil-flow visualization, however, showed separation 
occurring just ahead of the transition trip. Therefore, the height of the trip 
was ground off in progressive steps until no separation was apparent. The trip 
height was 0.025 inch when no separation occurred, and it was this height that 
was used for the results presented herein. A photograph of the leading-edge 
ring with the trip is shown in figure 2. 


Instrumentation 

Static-pressure orifices were located circumferentially on the cylinder, 
90° apart, at each of the 18-, 26-, and 54-inch stations (measured from the 
leading edge). There was also another orifice located 10 inches from the 
leading edge. (See fig. 1.) The static pressures were measured with a Statham 
pressure transducer (range from 0 to 1 psia) , mounted on an electrically actu- 
ated pressure scanning valve so that all the static orifices were read by the 
same transducer. This arrangement allowed a more accurate comparison to be 
made of the static pressures at various positions on the model. 

The pitot pressures obtained during the boundary-layer survey were meas- 
ured by three Statham pressure transducers manifolded together and connected to 
another electrically actuated pressure scanning valve. The three transducers 
had the following ranges: from 0 to 1 psia, from 0 to 5 psia, and from 0 to 

15 psia. This arrangement provided increased accuracy by permitting the pres- 
sure measurements to be taken from the gage which was nearest full-scale 
reading. A zero reading was taken on all pressure transducers before each run 
by using a reference pressure. 

Five swaged copper-constantan thermocouples were spot welded to the inside 
surface of the 0.12-inch-thick outside wall of the model at distances of 6, l4, 
22, 30, and 38 inches from the leading edge. The thermocouples were referenced 
to an insulated junction box outside the tunnel. The temperature of the insu- 
lated junction box was read on a 12-channel recording potentiometer. All the 
thermocouple and pres sure -transducer outputs were recorded on a high-speed ana- 
log to digital data recording system. 


Probes and Survey Apparatus 

Schematics of the boundary-layer total-pressure and total-temperature sur- 
vey probes used in the tests are shown in figures 3 and 4, respectively. The 
pressure probe had an outside diameter of 0.050 inch and was made of stainless- 
steel tubing flattened to a height of 0.007 inch with an opening height of 
0.003 inch at the tip. This tube was silver soldered into progressively larger 
tubing leading to the tubing (0.25- inch outside diameter) which formed the ver- 
tical shaft of the probe. 
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The temperature probe was made from a swaged thermocouple with a 0.0l4- 
inch outside diameter. The end of the swaged thermocouple was flared to approx- 
imately 0.020-inch diameter and two 0 . 007 -inch-diameter holes were drilled in 
the side wall for exits. The thermocouple was of 0.001-inch-diameter chromel- 
alumel wires, and the reference junction was placed in an ice hath outside the 
tunnel. 

The probe survey mechanism was driven by an electric motor and could posi- 
tion the survey probe in the boundary layer with an accuracy of ±0.0005 inch. 

The probe was electrically insulated from ground, and the position of the model 
surface was indicated by electrical contact of the probe with the model after 
equilibrium conditions were established in the tunnel. The survey mechanism 
also allowed the probe to be moved 10 inches in the streamwise direction while 
the tunnel was in operation so that various stations on the model could be 
surveyed. 


TESTS 


The tests were conducted at a free-stream Mach number of 6. The Reynolds 
number per foot was approximately 10 X 10^ obtained with a tunnel stagnation 
pressure of 525 psia and a stagnation temperature of 415° F ± 10° F. In order 
for the model wall to reach equilibrium temperature, the tunnel was operated 
for approximately 10 minutes before boundary-layer measurements were made. 
Total-pressure surveys of the boundary layer were made at seven different sta- 
tions (however, two surveys were made at the 40-inch station). These stations 
were along the top of the model at distances of 5> 6, 8, 11, 33* 37 > and 
40 inches from the leading edge. In addition to the total-pressure surveys, a 
total-temperature survey of the boundary layer was made at the 33-inch station. 
Surface temperatures and static pressures were measured during each survey. 
Schlieren photographs were taken in order to show the boundary-layer development 
for both natural and artificially produced transition. 

At the initiation of the test program, the model had to be alined with the 
stream flow. Figure 5(a) gives the static-pressure distribution around the 
model with the initial alinement. Also included in the figure are the angles 
of attack and yaw calculated from two-dimensional theory. Figure 5(*b) gives 
the same information after final alinement. The angle of attack calculated 
from the latter distribution was less than -0.05° and the angle of yaw was less 
than 0.10°. 


RESULTS AND DISCUSSION 


Turbulent Recovery Factors 

With the five thermocouples, wall temperatures were measured at various 
distances from the leading edge. The recovery factors were calculated from 
these measurements in order to check the validity of the assumption of zero 
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heat transfer. The recovery factors are defined by 


n 


T w - T& 




(1) 


and are plotted against x^ e in figure 6. The quantity T§ is the local 

static temperature at the edge of the boundary layer at any particular station. 
The average value of the turbulent recovery factor is seen to be approxi- 
mately 0.88. 


The radiation effect on the wall temperature was theoretically determined 
at the maximum radiation condition which occurred on the part of the model sur- 
face that can be seen through the schlieren windows. (See fig. 1.) Here the 
effect was found to decrease the wall temperature approximately 4° and decrease 
the recovery factor about 0.5 percent. Since this effect was not significant, 
no attempt to correct the wall temperatures was made. 

The plot in figure 6 shows a gradual decrease in recovery factor with dis- 
tance from the leading edge. This decrease is also shown in references 11 
and 12 for which the recovery factors were obtained at Mach numbers of 2.4 and 
from O .87 to 5*05, respectively. At the low Reynolds numbers for which a com- 
parison can be made, the recovery factors of the present test agree very well 
with those in reference 12 at a Mach number of 5«05* 


Total-Temperature Profile 

Total-temperature survey .- A total-temperature survey of the boundary 
layer as well as a pitot-pressure survey was made at the 33-inch station. The 
calibration of the temperature probe showing the variation of the probe recov- 
ery factor with Mach number and Reynolds number per foot was obtained from 
data taken in the Langley Unitary Plan wind tunnel at stagnation temperatures 
of 150° F and 230° F and is shown in figure 7- Also shown in the figure are 
the free- stream recovery factor obtained during the actual total-temperature 
survey and the interpolated curves necessary to go from the highest calibrated 
recovery factor to the free-stream recovery factor. The probe recovery factor 
is defined by 


T p - T 
n P = T t - T 


( 2 ) 


With this calibration and the Mach number and temperature -survey data at 
the 33-inch station (station 5} tables I and II), the total-temperature profile 
can be obtained by a quickly convergent iteration process of assuming a probe 
recovery factor rjp, finding the total temperature and Reynolds number from the 

assumed rjp, and then finding a new T]p with the Mach number and calculated 
Reynolds number. This process is repeated until the calculated T]p agrees 
with the assumed value. The final result of this iteration is shown in the 
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calibration plot (fig. 7) which gives the variation of Reynolds number with 
Mach number through the boundary layer. The total -temperature profile obtained 
in this way is shown in figure 8. As expected, several points in the boundary 
layer had total temperatures higher than the tunnel stagnation temperature. 

An indication of the accuracy of the temperature measurements can be 
obtained by finding the boundary-layer total-enthalpy thickness. Theoretically, 
the total-enthalpy thickness 0g should be zero for adiabatic wall conditions. 

The total-enthalpy thickness on a cylinder is defined as 



Figure 9 shows the integrand of equation (3) obtained from the pressure and 
temperature surveys at the 33-inch station and plotted against the height y. 
Good accuracy of the boundary-layer temperature probe is indicated inasmuch as 
the area under the curve, the total-enthalpy thickness 9 jj, is virtually zero. 

Comparison of measured temperatures with theory .- Since both a temperature 
and a pitot -pres sure survey were made at the 33-inch station, a comparison of 
the measured temperature-velocity relationship with various assumptions for this 
relationship can be made. One assumption, commonly used at low Mach numbers, 
is that the total temperature is constant through the boundary layer. Another 
assumption which was used in reference 13 is that the temperature distribution 
through the boundary layer is similar to the Crocco quadratic form (ref. l4) . 
This distribution is given by: 





where 



and 


D 


Tw 


(zero heat transfer) 


When equation (^) is transformed into its corresponding total-temperature dis- 
tribution, it is seen to be of the form: 
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(5) 




T 


t ,6 



+ b 


u 



+ c 


where 


a = 1 - 


aw 


t,5 


b = 


^aw T v 

t , 6 t , & 


= 0 


and 


c = 


L w 


L t,5 


(zero heat transfer) 


Another assumption that might be made for the temperature-velocity rela- 
tionship in a turbulent boundary layer is that it follows that of Crocco for 
the laminar boundary layer. This distribution is given by Van Driest in ref- 
erence 14 and in the nomenclature of this report is 


T t = 
t , 5 ^ 1 , 5 



JV 

T t,5. 



+ 





( 6 ) 


where 9^ and 0^ are coefficients which are given in reference l4 and 
depend only on Prandtl number and velocity ratio u/Ug. For use in this 
report a Prandtl number of 0.725 was assumed. 


A plot of boundary-layer total-temperature ratio as a function of the 
velocity ratio squared is given in figure 10. This figure shows the quadratic 
and Crocco laminar distributions and the assumption of a constant total temper- 
ature, all compared with the measured total temperatures. It appears from this 
plot that the quadratic temperature distribution would more adequately match 
the measured distribution than would the assumption of a constant total temper- 
ature or the Crocco laminar distribution. This agreement is especially true in 

the inner portions of the boundary layer ^i.e., for (u/Ug^ 2 < 0.7^. However, 

a comparison in figure 11 of the momentum profiles at the 55-inch station using 
the measured temperature distribution and the theoretical distributions shows 
that the assumption of a constant total temperature gives a momentum thickness 
very close to the value obtained from the measured temperature, whereas the 
quadratic and Crocco laminar distributions give values that are 15 percent 
greater and 16 percent less, respectively. 


A comparison of figures 10 and 11 indicates the necessity for a good 
theoretical turbulent temperature-velocity distribution to match closely the 
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temperature "bump" (that portion of the profile which exceeds the temperature 
T t,s) in tlie outer portions of the boundary layer (fig. 10). The constant 
total temperature does not correspond with either the inner or outer portions; 
but when the momentum thickness is determined, the errors are compensating and 
give a reasonably accurate value for the momentum thickness. The quadratic 
distribution, though corresponding to the temperatures in the inner portion of 
the layer, fails to account for the high temperatures of the bump and, as a 
result, gives a value of momentum thickness that is too high. On the other 
hand, the Crocco laminar distribution predicts temperatures that are consis- 
tently too high for all parts of the turbulent boundary layer and gives a 
momentum thickness which is too low when compared with the momentum thickness 
of the measured distribution. 


Reduction of the Pres sure -Survey Data 

In order to obtain the momentum thickness and velocity profile from the 
pitot-pressure survey, the temperature distribution must be known or a 
temperature-velocity relationship must be assumed. Since a total-temperature 
profile was measured only at the 35-inch station, it was necessary to assume a 
temperature-velocity relationship for the other survey stations. It was decided 
to use the temperature-velocity relationship obtained at the 33-inch station for 
the other stations also. It is realized that this relationship may be somewhat 
different at the other stations, but it seems certain that use of this particu- 
lar temperature -velocity relationship should give more accurate velocity pro- 
files than the commonly used assumption of a constant total temperature. As 
stated in the previous section, there was close agreement in the momentum thick- 
ness as determined by the two methods at the 33-inch station, but a comparison 
at the forward stations showed that they differ by as much as 5 percent. The 
method which gives the most accurate momentum thickness cannot be determined; 
however, in order to be consistent, the measured temperature-velocity relation- 
ship was also used to determine the momentum thickness. The measured- 
temperature-velocity relationship is shown in figure 10 in which the total- 
temperature ratio is plotted against the velocity ratio squared. 

An expression for the ratio of the velocity in the boundary layer to the 
velocity at the edge of the boundary layer is as follows: 


u M / l + °-3<S g Tt 

U S M SV1 + 0.2M 2 T t,5 


( 7 ) 


The Mach number distribution was obtained from the ratio of the pitot-survey 
pressure to the static pressure. The static pressure was assumed to be con- 
stant through the boundary layer. 

With the Mach number distribution and the Mach number at the edge of the 
boundary layer known, the other two unknowns of equation (7) (i.e., the ratios 
u/Ug and T^/T^ g) were determined by an iterative process. First, a total- 

temperature ratio T^/T^g was chosen and then a velocity ratio u/Ug was 
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calculated by using equation (7) • This calculated u/Ug was then used to find 
a new T^/T^g from figure 10. This process was repeated until agreement was 
reached between the two values of T^/T^g. The boundary-layer pressure, Mach 

number, and velocity ratios for each of the seven survey stations are presented 
in table I. 

With the total-temperature ratio T-g/T-^g and velocity ratio u/lJg known, 

the boundary-layer momentum thickness for the cylinder was calculated by using 
the following equation: 


9 = 



PU A _ u\ r + y 

P& U S V u &/ r 


ay 


( 8 ) 


where 


P _ ^~5 

P5 = t" 


The displacement thickness was calculated by using the fol- 


lowing equation: 


5* = 



_ dy 

P6 U 5/ r 


(9) 


The boundary-layer thickness 5p r was determined by plotting the boundary- 

layer survey pressures p' against the measured height y in the boundary 
layer. An example of one such plot is shown in figure 12. Points outside the 
boundary layer are at a nearly constant pressure. The intersection of a 
straight line drawn through these points and a straight line drawn through the 
last few data points just inside the edge of the boundary layer was taken to be 
the total thickness. The velocity ratio at the total thickness when found in 
this manner was approximately 0.998. 


Values of boundary- layer momentum and displacement thicknesses and of 
boundary-layer thickness based on pitot-pressure profiles are presented in 
table III, along with the skin-friction coefficients and Reynolds numbers. 


Turbulent Velocity Profiles 

The nondimen sional velocity profiles for the different stations are 
plotted in figure 13 as u/Ug against y/e. Power profiles, one-seventh and 
one-ninth, are also plotted to serve as a basis for comparison. The profiles 
indicate that n decreases somewhat with an increase in Reynolds number and 
becomes constant at the last three stations (fig. 13 (b)) where the profiles 
correspond to a one-eighth power. The profile at the 8-inch station 
(fig. 13 (a)) corresponds more closely to a one-ninth power which is a higher 
value of n for the measured value of Rq than that indicated in a correla- 
tion of n with Rq given in reference 15- It is unlikely that the larger 
value of n is associated with incomplete transition at the 8- inch station. 
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inasmuch as reference 1 6 indicates that natural transition should be completed 
at an approximate Reynolds number of 5 X 10^ or, in this case, 6 inches. 

In order to see whether the transition trip was causing any distortion in 
the profiles, a survey of the boundary layer was made at the 11-inch station 
with a smooth leading edge on the model. By a comparison in figure 13(c), the 
smooth-leading-edge profile is seen to correspond to some profile between the 
8- and 11-inch stations. This result is to be expected since transition 
occurred later on the smooth model, as shown by the schlieren photographs 
(figs. l4(a) and (b))j therefore, the effective Reynolds number of this profile 
would be lower than that of the 11-inch station with a transition trip. Since 
the profiles with a transition trip have the same general shape as the profile 
with natural transition, it is believed that there was no distortion effect due 
to the transition trip. 


It is of interest to note here that even though the profiles correspond to 
a changing power profile with Reynolds number, the shape factor 5*/0 was 
found to be relatively constant except for the 5-inch station which is thought 
to be in the region of transitional flow. These shape factors for the dif- 
ferent stations are presented in table III. An average value for 8*/0 was 
found to be 16.9 with a variation of less than 2 percent. This nearly constant 
shape factor is to be expected since, for a power-law profile of the form 
u / v 

— = (c I j it is relatively independent of the value of n. This relative 

U5 \o/ 


independence is shown in figure 15 in which the shape factor 6*/0 is plotted 
against n. Little variation of 5*/0 with n is shown, especially at the 
higher values of n. The transverse curvature of the boundary layer is indi- 
cated by the factor 5/r. The effect of this factor on the shape factor 5*/0 
is also shown in figure 15* The effect of 5/r is seen to be negligible for 
the maximum 5/r encountered in this investigation, that is, 5/r « 1/5. 


Theoretical values of 5/0 are shown in figure 15 • They are seen to vary 
considerably with both 5/r and n. The method from which the theoretical 
curves of figure 15 were calculated is presented in the appendix. 


Skin Friction 

One of the main objectives of this investigation was to determine average 
skin-friction coefficients for the turbulent boundary layer by the momentum 
method at a Mach number of 6 and adiabatic-wall conditions. 

Longitudinal Mach number distribution .- The Mach number distribution along 
the top of the model obtained from probe total pressures, static pressures, and 
the tunnel stagnation pressure is shown in figure 16. The Mach number increases 
from about 5-95 at 5 inches to 6.02 at 40 inches. The skin-friction coeffi- 
cients should be obtained when the flow conditions along the edge of the 
boundary layer are constant . In reference 6 Wilson compensated for a change in 
Mach number from the front of the model to the survey station by integrating his 
momentum thicknesses across a fictitious plane removed from the survey plane. 

The conditions outside the boundary layer of this fictitious plane correspond 
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to the average Mach number along the entire model, that is, Mg. In the nomen- 
clature of this report, Wilson's equation for the momentum thickness at the 
fictitious plane in terms of the conditions at the survey plane and the Mach 
number outside the boundary layer at the two planes is 


7+1 



Instead of tediously integrating 9 from the survey data, a procedure 
similar to that of Wilson's was followed; the ratio 9/9 for different values 
of Mg was calculated by assuming the following velocity distribution: 



The momentum-thicknessratio is plotted against the ratio of the Mach numbers 
-of the two planes Mg/Mg in figure 17- The algebraic average Mach number over 
the entire cylinder was calculated to be 5-98 (i.e., Mg = 5*9' 8j. Since the 

Mach number along the model did not vary greatly, this correction to 9 was 
very small, being less than 1-5 percent. 

Determination of virtual origin .- Since all theoretical analyses are based 
on the assumption of fully developed turbulent flow from the leading edge, the 
region of laminar flow in this test has to be taken into account and, as 
explained in reference 2, can be done by finding the fictitious (virtual) origin 
of the turbulent boundary layer. The Reynolds number based on the length from 
the virtual origin is the effective Reynolds number of the turbulent boundary 
layer . 

The various methods for determining the position of this origin are 
described in reference 2. Two of these methods were applied to the present 
data. The first method used was the extrapolation of the momentum thickness, 
measured in the turbulent portions of the boundary layer, to zero. This extrap- 
olation was done by assuming a growth of the boundary layer of the following 
form: 


R 0 - K( R le - R le-x) C 


( 12 ) 


This mathematical curve was fitted to the experimental data by the method of 
least squares. Values obtained for R^ e _ x an ^ C were I .56 x 10° and 0.679* 

respectively, and when substituted into equation (12) yield a virtual origin 
that is about I .87 inches downstream from the leading edge. 
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The second method used was that of Rube sin, Maydew, and Varga (ref. 17)* 

In this method the same type of power equation is assumed for the growth of the 
boundary layer, but the constant C is obtained by averaging the values pre- 
dicted by several theories. For the conditions of this investigation, the aver- 
age value of C used was 0.800, the data were plotted on log-log paper, and 
the value of R j e _ x necessary to match the slope of the line C = 0.800 was 

obtained. This procedure yielded a value of R-^e-x “0*0^ x 10^* There- 
fore, the virtual origin by this method is 0.05 inch upstream from the leading 
edge. 


A plot of Rq against the Reynolds number based on the distance from the 
leading edge R^ e is shown in figure 18(a) . Station 1 (5-inch station) is 

thought to be in the region of transition (note data symbol labeled "transi- 
tional flow" in fig. 18(a)). Besides the fact that this data symbol does not 
fall in line with the data from the other stations , the schlieren photo of fig- 
ure l4(b) indicates that this station is in the transitional flow region. Con- 
sequently, the momentum thickness from this station was not used in the deter- 
mination of the virtual origin of the turbulent boundary layer. 

Figure 18(b) shows a plot of Rg against the Reynolds number based on the 

distance from the virtual origin as found by the two methods. Since it is not 
known which of the two methods is more accurate, it is considered that the 
actual virtual origin could be at any point between the calculated virtual 
origins. 

Average skin-friction coefficients .- The average skin-friction coefficients 
can be determined from the equation 


c p = 2 ?£ (15) 

P R x 

if the virtual origin is known so that the effective Reynolds number R x can 

be determined. A plot of the average skin-friction coefficients at various 
Reynolds numbers, determined by using the calculated virtual origins, is shown 
in figure 19. Also shown in the figure are four turbulent skin-friction 
theories. The variation of the average skin-friction coefficients due to the 
difference in the virtual origins is seen to be rather large at the front sta- 
tions (20-percent variation from the average at the most forward station pre- 
sented). This variation decreases with increasing Reynolds number and at the 
back stations is only about 3 percent. 

Since the present experimental values of Cp are for a cylinder, the 
theoretical values of Cp for a flat plate should be modified to take into 
account the effect of the lateral surface curvature. An increase in Cp of 

about 1 percent because of surface curvature is shown in reference 18 for a 
ratio of boundary-layer thickness to cylinder radius of 0.l8. This ratio of 
0.18 is the largest ratio of boundary-layer thickness to cylinder radius 


occurring in this investigation. The correction at any forward station would 
be’ less. For the incompressible case, references 19 and 20 also give a cor- 
rection of about 1 percent at the last station. Since the effect of curvature 
is seen to be so small, no attempt was made to correct the skin-friction 
coefficients. 

The average skin-friction data are replotted in figure 20 as ratios to the 
incompressible skin-friction coefficients. The incompressible skin-friction 
coefficient was obtained from the well-known Karman-Schoenherr equation for 
average skin-friction coefficient in incompressible flow. From figures 19 
and 20 it can be seen that it would be difficult to say that the experimental 
average skin-friction coefficients correspond to any one particular theory 
because of the uncertainty at low Reynolds numbers. At the high Reynolds num- 
bers where the exact location of the virtual origin is not as important, the 
coefficients are between the values predicted by the theory of Wilson and 
Van Driest and the T' method of Sommer and Short. 

Along with the present data, other authors' experimental values of average 
skin-friction coefficients (refs. 1, 6, 17 ; 21, and 22) are presented in fig- 
ure 21. Here the ratio of the skin-friction coefficient to the Karman- 
Schoenherr incompressible coefficient is plotted against Mach number at a con- 
stant Reynolds number of 10 X 10^. Since the choice of a Reynolds number of 

10 x 10^ for figure 21 was an arbitrary one, it should be noted that the rela- 
tion between the present data and the theories would change for any other 
Reynolds number chosen. This changing relation can be seen by referring to 
figure 20. 

As figure 21 shows, the average skin-friction data available at this time 
were obtained at Mach numbers less than 4 and, therefore, comparison with the 
present data is difficult. Local skin-friction data, however, have been 
obtained at higher Mach numbers by Korkegi (ref. 23) and Matting, Chapman, 
Nyholm, and Thomas (ref. 24). Therefore, local skin-friction coefficients were 
obtained from the present data for comparison. 

Local skin-friction coefficients .- The local skin-friction coefficient for 
the present data was obtained from the rate of growth of the measured momentum 
thickness by differentiating the following equation of the curves in fig- 
ure 18 (b) which are of the form 


Rg = KR X C (l4) 

and then using the equation 

dRs 

c f - 2 (15) 

In figure 22, local skin-friction data are plotted against Mach number in 
a manner similar to that of figure 21 for the average skin-friction data. 
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Experimental local skin-friction coefficients obtained by other authors were 
also plotted in figure 22 and these data can be found in references 17 and 
23 to 26. Korkegi obtained local skin-friction coefficients from force-balance 
measurements at a Mach number of 5-8 and Reynolds numbers from 2-3 X 10^ to 
3.9 X 10^. The data were extrapolated, as is to be explained, to R x = 10 x 10^ 
for figure 22 of this report. Matting and coauthors (ref. 24) also used the 
skin-friction balance to obtain local coefficients at Mach numbers of 2.95 
and 4.2 in air. Then, in helium flow, local skin-friction data were obtained 
at equivalent air Mach numbers of 4.2, 6.7, and 9-9- The Reynolds number var- 
iation was from 2 x 10^ to 100 x 10^. It was found that the data from the 
measurements in air and helium agreed very well at the Mach number of 4.2, but 
no comparison is available at the higher Mach numbers. 

Figure 22 shows how the present data compare with these two sets of data. 
Korkegi' s coefficient at a Mach number of 5-8 is about 12 percent higher than 
the present data, whereas the level of Matting's high Mach number data obtained 
in helium is considerably lower than that of the present data. The plotted 
point for Korkegi 's data was obtained by extrapolating the experimental values 

of Cp and the effective Reynolds numbers (R x ,max ~ 5*5 X 10^), as determined 

in reference 2, to the Reynolds number of 10 x 1C)6. (Other extrapolated points 
in figs. 21 and 22 were obtained in a similar manner.) Korkegi 's data are such 
that considerable variation of the extrapolation is possible, and the bar 
through the plotted point gives an indication of this variation. 

In reference 2, a comparison of data available at that time indicates that 
the Sommer and Short T' method appeared to provide the most accurate estimate 
of skin friction within the ranges of Mach number, Reynolds number, and heat 
transfer of the available experiments. The present data, however, casts doubt 
as to which theory will provide reliable estimates of skin friction at Mach num- 
bers greater than 4. Additional experimental research on turbulent skin fric- 
tion is required at hypersonic speeds. 


CONCLUSIONS 


An investigation of the turbulent boundary layer on a hollow cylinder at a 

Mach number of 6, at Reynolds numbers up to 33 X 10^, and for adiabatic-wall 
conditions leads to the following conclusions: 

1. For the turbulent boundary layer at a Mach number of 6, the temperature 
recovery factor was found to be approximately 0.88. This value is in very good 
agreement with values of other experiments. A gradual decrease in recovery 
factor with an increase in Reynolds number was obtained which is also in accord 
with findings of other investigations. 

2. The experimentally determined total-temperature-velocity relationship 

did not agree with any of the three theoretical distributions considered: the 

assumption of a constant total temperature, the quadratic distribution used by 
Donaldson in NACA RM L52H04, and the Crocco laminar distribution used by 
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Van Driest in NACA TN 2597* However, the momentum thickness calculated hy 
assuming a constant total temperature agreed well with the momentum thickness 
determined with the experimentally measured temperature distribution. The 
absence of a good theory for the temperature-velocity relationship in turbulent 
boundary layers is a problem which merits more attention. 

3. A comparison of measured velocity profiles with a power-law velocity 

profile of the form = 

U S 

showed a decrease in the exponent n from 9 at a Reynolds number of about 

6 X 10^ (8-inch station) to 8 at a Reynolds number of approximately 32 X 10^ 
(40-inch station). The value of n at the lower Reynolds number was somewhat 
greater than that indicated by the correlation of experimental data in NAVORD 
Rept. 3854 by Persh. The boundary-layer shape factor, however, was found to be 
relatively constant at a value of approximately l6*9; which agrees well with 
the theoretical shape factor for a power -velocity profile . 

4. There was considerable uncertainty in the average skin-friction coef- 
ficient at low Reynolds numbers as a result of the difference in the position 
of the virtual origin as calculated by the methods of least squares and of 
Rubesin, Maydew, and Varga (NACA TN 2305)* This uncertainty made it difficult 
to compare the data with theories in the low Reynolds number range. At the 
higher Reynolds numbers, however, the uncertainty in the value of the skin- 
friction coefficient was less and the coefficients were between the values pre- 
dicted by the theory of Wilson and Van Driest and the T f method of Sommer and 
Short. Exact establishment of the hypersonic turbulent skin-friction level 
requires additional experimental research. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 19* 1965* 


(where 


is the bound ary -layer thickness) 
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APPENDIX 


EQUATIONS FOR DETERMINATION OF 5*/e AND 5/e FOR 
POWER-LAW BOUNDARY-LAYER VELOCITY 
PROFILES ON A CYLINDER 


The nondimen sional equations for 6* and 0 for the turbulent boundary 
layer on a cylinder are as follows: 


T • 1 (‘ • sSfe)(‘ * H)- i 

i ' X‘ S(‘ - «(' * 


(Al) 


(A2) 


If the assumptions of a constant total temperature and pressure through the 
boundary layer and a power-profile variation of velocity through the boundary 

■q / y Y^* / ^ 

layer of the form - — = [ — ] are made, then exact solutions for equa- 

u 6 \ 5 / 

tions (Al) and (A2) can be obtained in the following way. For a constant total 
temperature through the boundary layer, the following energy equation is true 
for all points in the boundary layer: 


U 1 


+ c p T - c p T t “ c p T t,6 


The density distribution is 


P6 _ T t,6 / T t,5 A / u \ 2 
P " t 5 ' t b 




By setting z = — = 


becomes 


U* 


A) 


,1/n 


a = 


Tr Tg 

and b = 1 - , equation (Al) 


L t,5 


T. 


t,5 


$ ■ £ (‘ - nt # * ? -)-' 1 


dz 
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APPENDIX 


When the integral is split and simplified, this equation becomes 


5 


1 


na 


f- 

J 0 1 


z n , , 15 S „ 0 

dz + — na 

2 r r 


bz2 


I 1 - 

J 0 1 - 


2n 
bz2 


dz 


(A3) 


The first two terms are those obtained for a flat plate in reference 27. 


If n takes on integer values, the exact solutions to the integral equa- 
tion (A3) for 5*/5 can be written in series form. The series solution for 
values of n which are odd integers can be written in the following summation 
form: 


6 * 

5 


1 


+ 


i® + 

2 r 


n a 
2 t (n+l)/2 



b 

T + l0g e 


+ 




2i + 1 



1 + /b 
1 - /b 


w 


In finding a similar equation for the momentum thickness, the quantities 
a, b, and z are substituted into equation (A2) and the following equation 
is obtained: 


0 

5 


na 



(A5) 


When the values of n are restricted to odd integers, the results of the inte- 
grals in equation (A5) can be written in summation form as follows: 


0 

5 


n 


2 (n+l)/2 


n-1 

2 


2 - log e a + 


1 i 1 - /b V b^~ 
ft e 1 + ft ' 1(21 + 1 ) 


, 6 n a 
r 2 b n 


n-1 


1 , i n _ 1 - {b b^ 

^ l°ge a - - °Se x + ^ - 2, (1 + l)(2i ♦ 


1) 


(a6) 


Equations (a 4) and (a 6) were used to obtain the values of 5*/9 and 5/0 in 
figure 15. 
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TABLE I.- BOUNDARY- LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS 


(a) Station 1: x Je = 5-0 in.; T 5 = 109° R; T t = 878 ° R; Mg = 5- 

P 5 = 0-353 psia; and p t = 523 psia 


y> 

in. 

P5 

p' 

mm 

U 

u 5 

■ 

0.2788 

0.259 

0.5855 


.2102 

.305 

.6525 


.1618 

.353 

• 713 


.1192 

.4i6 

.778 

.0121 

.0951 

.469 

.824 

.0162 

.o 8 o 4 

• 512 

.855 

.0218 

.0641 

.576 

.894 

.0288 

• 0571 

.611 

.9115 

.0372 

.0522 

. 64 o 

.9345 

.0520 

.0415 

.720 

• 9535 

. 0668 

• 0335 

.803 

.974 

•0793 

.0294 

• 859 

• 984 

.096^ 

.0245 

• 94i 

.994 

.1172 

.0219 


•9995 

.1387 

.0212 

1.013 

1.001 

.1637 

.0210 

1.017 

1.0015 


.0209 

1.019 

1.002 























TABLE I.- BOUNDARY -LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 

Ob) Station 2 : x Ze = 6.0 in.; Tg = 108° R; T t = 872 ° R; Mg = 5 - 

p 6 = 0.349 psia; and p^. = 523 psia 


y, 

in. 

P5 

p' 

M 

Mg 

U 

U 5 

0.0035 

0.2404 

0.282 

0.620 

.0045 

•2354 

.286 

.627 

.0086 

.2079 

.307 

•655 

.0154 

.1328 

• 393 

.756 

.0237 

.0929 

.474 

•8275 

.0365 

• 0642 

• 575 

.893 

.0490 

.0556 

.619 

•9155 

.0617 

.0476 

.671 

•9365 

.0870 

•0353 

.782 

•970 

.1122 

.0283 

.874 

.986 

.1375 

.0235 

•959 

.996 

1 -1501 

.0226 

.980 

•998 

.1711 

.0219 

•996 

•999 

.2027 

.0217 

1.000 

1.000 



TABLE I.- BOUNDARY- LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 


(c) Station 3: x le = 8.0 in.; Tg = 108° R; T t = 871° R; M g 

p & = 0.3^ psia; and p t = 519 psia 



1.0005 











TABLE I.- BOUNDARY-LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 


(e) Station 4: 


x le = 11.0 in.; Tg = 108° R; T t = 872 ° R; Mg = 5- 
Pg = 0-345 psia; and p^ = 523 psia 


y, 

in. 

£5 

P' 

M 

Mg 


0.0035 

0 . 253 ^ 

0.273 

0.607 

• 0035 

.2480 

.277 

.613 

.0068 

.2349 

. 286 

.626 

.0112 

.2051 

• 309 

.658 

• 0155 

.1501 

.367 

.7285 

.0236 

.1221 

. 4 io 

.772 

• 04 o 5 

.0935 

.472 

.826 

.0507 

.0849 

.496 

.8445 

.0616 

.0761 

• 525 

.864 

.0869 

.0598 

• 595 

• 9035 

.1123 

.0473 

.671 

.9365 

•1573 

.0381 

• 750 

.962 

.1645 

.0311 

.831 

.980 

.1906 

.0267 

.899 

.9895 

.2198 

.0239 

• 951 

• 9955 

.2385 

.0227 

.976 

.998 

.2697 

.0219 

•995 

• 999 

•2935 

.0219 

.994 

•9995 




TABLE I.- BOUNDARY-LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 


(f) Station 5: x Ze = 33.0 in.; Tg = 106° R; T t = 871° R; Mg 

Pg = 0.325 psiaj and p^. = 523 psia 


y > 

in. 

P6 

p' 

M 

Mg 

U 

0.0035 

0.2484 

0.274 

0.608 

• 0035 

.2478 

.274 

.608 

.0065 

.2405 

.279 

.6155 

.0085 

.2314 

.285 

.625 

.0166 

.1845 

• 324 

•6775 

.02^9 

.1568 

.354 

.714 

•0331 

.1442 

• 371 

• 733 

.0444 

• 1337 

.386 

.749 

•0570 

.1261 

•399 

.761 

.0697 

.1194 

.411 

.77^ 

.0950 

.1049 

.440 

.800 

.1202 

.0936 

.467 

.822 

.1V73 

.0828 

.498 

.846 

•1853 

.0701 

.543 

• 875 

.2255 

.0590 

• 593 

.9025 

.2803 

.0468 

.669 

• 9355 

.3309 

.0381 

.743 

.960 

.3815 

.0312 

.822 

.978 

• 4375 

.0258 

.905 

.990 

.4824 

.0229 

.961 

.9965 

• 5331 

.0217 

.988 

•999 

.5964 

.0214 

•995 

• 9995 

.6615 

.0213 

•997 

1.0000 



TABLE I.- BOUNDARY-LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 


(g) Station 6: * le = 57.0 in.; T 6 = 105° Rj T t = 867° R; Mg = 6.02 

Pg = 0.324 psia; and p^ = 523 psia 


y> 

in. 

P5 

P f 

M 

Mg 

U 

u 5 

0.0035 

0 . 2711 

0.260 

0.586 

•0035 

•2731 

• 259 

.586 

•0053 

.2637 

.264 

• 593 

.0094 

.2434 

.277 

.613 

.0166 

.1946 

• 315 

.666 

.0223 

.1620 

.348 

•707 

.0293 

.1481 

.366 

.728 

.0434 

.1340 

.386 

•749 

•0579 

.1247 

.401 

• 764 

•0715 

.1169 

.415 

•777 

.0968 

• 1037 

.442 

.8005 

.1222 

.0920 

.471 

.826 

.1474 

.0821 

.500 

.847 

. 2009 

.0646 

. 566 

.888 

.2499 

.0528 

.628 

• 9195 

•3013 

.0430 

.697 

.946 

.3626 

.0342 

• 784 

•970 

.4277 

.0276 

.874 

.986 

.4931 

.0232 

•955 

•996 

• 5522 

.0215 

•993 

•999 

.6178 

.0212 

1.000 

1.000 

.6813 

.0212 

•999 

1.000 







TABLE I.- BOUNDARY-LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Continued 

(h) Station 7: x Ze = 40.0 in.; T & = 106° R; T t = 871° R; Mg = 6. 

p & = O.326 psia; and p^ = 524 psia 


y, 

in. 

P5 

P' 

M 

Mg 

u 

U 6 

0.0035 

0.2464 

0.276 

0.612 

•0037 

.2454 

.276 

.612 

.0066 

.2368 

.282 

.620 

.0097 

.2229 

.292 

.6345 

.0163 

.1830 

.326 

.680 

.0221 

.1609 

• 350 

.709 

.0288 

.1467 

.368 

•730 

.0432 

• 1334 

.387 

• 750 

.0563 

• 1255 

.4oi 

.7635 

.0684 

.1188 

.412 

•774 

.0948 

.1058 

.438 

•7975 

.1188 

.0950 

.464 

.820 

.1438 

.0863 

.488 

.839 

.1844 

.0734 

•531 

.868 

.2230 

.0627 

.576 

.894 

.2707 

.0523 

.633 

.9215 

.3230 

.0434 

.696 

• 946 

.3864 

.0346 

.780 

.9695 

.4497 

.0281 

.867 

.985 

• 5131 

.0236 

.949 

•995 

.5762 

.0216 

.991 

•999 

.6015 

.0213 

.998 

1.000 
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TABLE I.- BOUNDARY-LAYER PRESSURE, MACH NUMBER, AND 
VELOCITY DISTRIBUTIONS - Concluded 

(i) Station 7: x le = 40.0 in.; T & = 104° R; T t = 86l° R; Mg 

p & = 0.525 psia; and p t = 525 psia 


Y , 

in. 

P& 

P' 

M 

“8 

u 

U 5 

0.0055 

0.2469 

0.275 

0.610 

.0054 

.2475 

• 275 

.610 

• 0055 

.2475 

.275 

.610 

.0085 

.2550 

.284 

.625 

.0107 

.2152 

.299 

.644 

.0178 

• 1757 

• 355 

.6895 

.0222 

• 1591 

• 352 

.711 

.0504 

.1452 

• 373 

•735 

.0446 

.1500 

• 393 

.756 

• 0599 

.1208 

.408 

• 7705 

.0870 

.1071 

.455 

•795 

.1102 

• 0975 

.458 

.8145 

.1496 

.0828 

.498 

.846 

.1880 

.0706 

• 541 

.874 

.2492 

.0555 

.615 

.9125 

• 5159 

.0454 

.695 

.945 

.5925 

.0551 

.798 

• 973 

.4555 

.0270 

.886 

.988 

.5168 

.0228 

.963 

.9965 

.5852 

.0211 

1.002 

1.0000 

.6529 

.0210 

1.006 

1.0000 

.6979 

.0210 

1.004 

1.0000 



TABUS II.- TEMPERATURE SURVEY 


^Station 5 (test 20 ) : 


= 33-0 in. ; Mg 


6.02; and 



y , 

in. 

^t,5 

(a) 

T P 

T t 

0.0100 

860 

773 

801 

.0100 

860 

775 

803 

.0143 

857 

775 

8 o 4 

•0159 

856 

776 

806 

.0184 

857 

778 

807 

.0224 

861 

782 

812 

.0310 

866 

784 

815 

• 0372 

855 

778 

809 

.0435 

860 

782 

8 l 4 

.0561 

860 

785 

816 

.0687 

851 

781 

813 

• 1235 

855 

796 

830 

.1445 

867 

815 

852 

.1961 

865 

819 

856 

.2496 

858 

822 

860 

.2983 

859 

829 

867 

• 3466 

860 

837 

872 

. 4 i 4 o 

864 

846 

872 

•4753 

859 

844 

862 

■ 5403 

859 

847 

859 


a Measured in settling chamber; assumed to be T^. g. 
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TABLE III.- REDUCED DATA FROM PRESSURE -SURVEY MEASUREMENTS 


to 
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ON b O OJ OJ 

OJ to to to 

VD o 

rH 
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ITN H 
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H 
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t— CO -4 

VO ON ON 
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^ 3 3 £ ^7 

Os ^ 
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1 

3 
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diameter of 0.075 inch. 





Figure 3*- Schematic of boundary-layer pressure probe. 
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velocity ratio squared. 




Figure 11.- Momentum profiles by various boundary-layer temperature distributions 

at 33-inch station. 
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(c) 8- and 11-inch stations 


including 11-inch station for smooth leading edge 
Figure 13.- Concluded. 




(a) Natural transition. L-65-71 

Figure lk.- Schlieren photographs of boundary -layer development. 
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(b) Artificial transition. 
Figure l4.- Concluded. 


L-65-72 
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Variation of ratio of average skin friction to incompressible skin friction with 
Reynolds number. Mg = 5*9^- 
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Rubesin and coauthors (ref. 17? Rj = 6 x 10°) 
O'Donnell (ref. 21? =6* 10 6 ) 

Chapman and Kester (ref. 1) 

Brinich and Diaconis (ref. 22) 

V/ilson (ref. 6) 

Present data, virtual origin by least squares method 
Present data, virtual origin by method of reference 17 
Solid symbols extrapolated from R x ^ max 
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Figure 21.- Variation of various theoretical and experimental average skin-friction coefficients 

with Mach number. L = 10 x 10^. 
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♦ Korkegi (ref. 23? » 5.5 * 10 ) 

A Matting and coauthors (ref. 24) 

Coles (ref. 25 J R^. mx = 8 x 10^; 

Q Shutts and coauthors (ref. 26) 

O Present data, virtual origin by least squares method 

□ Present data., virtual origin by method of reference 17 

Solid symbols extrapolated from R^. max 

Flagged symbols indicate data obtained in helium 
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